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We report on the status of the Monte Carlo event generator WOPPER. Version 1.1 of WOPPER describes four 
fermion production at LEP-II and beyond with leading logarithmic radiative corrections in the double pole 
approximation. These approximations are appropriate for almost all practical purposes, but the inclusion of these 
finite width effects and radiative corrections is nevertheless indispensable for LEP-II physics. 
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1. Introduction 

The Monte Carlo event generator WOPPER [1] for 
four fermion production through resonances 
at LEP-II and beyond is the latest addition to the 
Darmstadt/Siegen family of Monte Carlo event 
generators [2]. It is a true event generator that 
generates a sample of unweighted events which 
can be used directly in detector simulations for 
experiments at LEP-II and future linear e'^e~ col- 
liders [3] . The distinguishing features of WOPPER 
are: off-shell pair production in the double pole 
approximation and resummation of the leading 
logarithmic (initial state) radiative corrections. 

2. Requirements 

While we still need to add a direct observation 
of the triple gauge vertices (TGV) to the over- 
whelming indirect evidence for their existence, 
any potential deviation from the standard model 
values will in all likelihood be very small [4]. 
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Thus an observation of these small anomalous 
couplings in the permille range will only be pos- 
sible at an high luminosity collider beyond LEP- 
II, if all standard model corrections are know to 
an even better precision. This is only possible 
if several independent semi-analytical and Monte 
Carlo programs are available, which include (and 
agree on) all important contributions. 

In the nearer future, the prime physics objec- 
tive of LEP-II will be a precise measurement of 
the mass, which will provide an important 
cross check of the standard model and help to 
constrain possible physics beyond. For example, 
reducing the error below lOOMeV could close the 
window for light Higgses (and therefore the min- 
imal supersymmetric standard model) using con- 
straints from the electroweak radiative correction 
parameter Ar, as reported at this conference [5]. 
Proper accounting for finite width effects and ra- 
diative corrections is of crucial importance for this 
measurement. This calls for reliable calculations 
and Monte Carlo event generators. 

3. Features of WOPPER 

3.1. Radiative Corrections 

WOPPER concentrates on the gauge invariant 
subset of radiative corrections which is phc- 
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nomenologically most important. The leading 
logarithmic (a" ln"(s/TOg)) radiative corrections 
from the initial state leptons are resummcd to all 
orders, including the exponentiation of soft pho- 
tons. 

Because of the t-channel i^-exchange diagram, 
it is not possible to unambiguously separate the 
complete 0{a) initial state radiative corrections 
in a gauge invariant way. The restriction to 
the leading logarithmic contribution is therefore 
phenomenologically reasonable and theoretically 
sound. 

While the resummation of terms beyond C(q;^) 
is numerically irrelevant for the projected LEP-II 
statistics, the fully resummed form can be ob- 
tained at no extra cost in a Monte Carlo event 
generator and has the desirable side effect of an 
unproblematic probabilistic interpretation. All 
finite order approximations suffer from the so 
called ko problem: for some infrared cut ko on 
the emitted photon energy, the cross section 

a«fl-0(l).-ln4ln^Vao (1) 

in the channel with no emission of a photon be- 
comes negative. In the resummed form, however, 
the big negative term in (1) exponentiates and 
leads to a small but positive cross section. There- 
fore the generated event samples are physically 
meaningful and the corresponding cross sections 
do not depend on the soft photon cut introduced 
by the experimental resolution. 

Besides the theoretical advantage of keeping 
only the gauge invariant leading logs, these cor- 
rections have the benefit of leading to a factorized 
cross section. Denoting the phase space variables 
under consideration collectively by [PS], the cross 
section reads 
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where the [PS'] depend on [PS], and the energy 
fractions a;^ of the initial state leptons. For ex- 
ample, we have s' = x~^x~s for the center of mass 
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Figure 1. Effect of radiative corrections on the 
total cross section near threshold. 



energy squared. This cross section can be imple- 
mented easily in an event generator. Here the 
electron distribution functions D{x;iJ.^) obey the 
DGLAP renormalization group equation: 
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Prom the cross section (2) it can be seen that 
the important effects are either universal or of a 
simple kinematical origin: the shape of the total 
cross section will be shifted towards higher en- 
ergies because of the energy loss from radiated 
photons ("radiative tail", cf. figure 1). Parti- 
cles which would be back-to-back due to momen- 
tum conservation in non-radiative events will be 
acollinear in radiative events. 

The Monte Carlo implementation of (2) in 
WOPPER first solves a suitably infrared regular- 
ized version of (3) with a universal photon shower 
Monte Carlo and use the sample corresponding to 
this solution to define an effective center of mass 
system (CMS) for each event after radiation of 
initial state photons. In this new CMS a rela- 
tively simple Born typo e+e^ 4/ event can be 
generated, which will then be boosted back to the 
laboratory frame. This final boost incorporates 
all kinematical effects, like acoUinearities for the 
intermediate Ws, etc. 

While this approach is more than adequate for 
all inclusive distributions, there remain two areas 
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Figure 2. E'^ spectrum at ^/s = 170 GeV, com- 
paring WOPPER's LLA prediction with wwf 's 0{a) 

prediction. An angular cut for photons around 
any charged particle in the final state has been 
applied. 
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Figure 3. Same as figure 2, but for the photonic 
transversal momentum p^. 



where further progress is needed: large photonic 
and final state radiation. 

By its very definition, the leading logarithmic 
or pole approximation is applicable to coUinear 
radiation and inclusive spectra because these are 
dominated by collinear radiation (cf. fig. 2, where 
the LLA is compared to a 0{a) calculation). But 
for large photonic the LLA is essentially an un- 
controlled approximation. However, wc can use 
existing 0{a) Monte Carlos [6] to gauge the nu- 
merical accuracy of the LLA even for large pj. It 
turns out that the pole approximation reproduces 
the full 0{a) calculation surprisingly well in the 
projected LEP-II energy regime (cf. fig. 3). At 
500GcV on the other hand, the large region is 
overestimated (cf. fig. 4). 

The other problematic issue concerns final state 
radiation. While the figures 2, 3 and 4 confirm 
that a reasonable angular cut around any charged 
particle in the final state will remove all effects of 
final state radiation, are more complete Monte 
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Figure 4. Same as figure 3, but at higher energy 
\/i = 500GeV. 
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Carlo treatment is of course desirable. 

3.2. Finite Width 

For detailed background studies, it is of course 
necessary to study all non-resonant diagrams con- 
tributing to e+e~ Af at LEP-II energies [7]. 
For the identification of W'^W~ pairs, experi- 
mental invariant mass cuts will have to be ap- 
plied, however. After such cuts, the contribution 
of non-resonant diagrams goes down rapidly [7] . 

For the sake of efficiency, version 1 . 1 of WOPPER 
implements the four fermion final states therefore 
in double pole approximation. This is equivalent 
to keeping the two so-called "signal" diagrams. 

For all practical purposes, the W^^fi decay into 
light quarks only (6's are still relatively light and 
because of the small \Vcb\'^ rare). Thus the un- 
physical polarizations decouple and the cross sec- 
tion factorizes 
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into "oflE-shell" production aoS-shen{s; s+, S- 
Breit-Wigner propagators 
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and decay widths Tw{s±)- This factorizcd form 
is again very convenient for implementation in 
a Monte Carlo event generator. The resulting 
cross section, displaying the typical smearing of 
the threshold, is depicted in figure 5. 

3.3. Parton Showers, Fragmentation and 
Hadronization 

The description of scmileptonic and hadronic 
W^W^ events is of course incomplete without 
proper accounting for the hadronization of the 
quarks in the final state. Working interfaces of 
WOPPER to both major final state parton shower 
and fragmentation models JETSET 7.4 [8] and 
HERWIG 5.5 [9] are implemented as of version 1.1. 
Therefore the events generated by WOPPER can 
be fed immediately to detector simulation Monte 
Carlos. 
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Figure 5. Effect of the finite W"^ width on the 
total cross section near threshold. 



4. Sample Application 

At high energies, the trans versally polarized 
W^'s from the f-channel diagram are dominat- 
ing the cross section. From a physics perspective 
however, the longitudinally polarized T4^^'s and 
the contribution from the s-channel diagram are 
both far more interesting. The former because 
the longitudinal M^^'s are a direct manifestation 
of electroweak symmetry breaking and the latter 
is where the three gauge boson vertex is to be 
measured. 

It is therefore important to be able to recon- 
struct a sample of longitudinally polarized T^^'s. 
Fortunately, the V — A decay of ff' is self 

analyzing, if we can measure the angular decay 
distribution: 

a+,o,-= /dcosr P+,o,-(cosr)-^ (6) 

Here the P+fl,- are simple polynomials [10]. 
In the case of on-shell W^^s without radiative 

corrections, the decay angle can be reconstructed 
in semileptonic decays by measuring the energy 
of the charged lepton only [10] 

2Ei — Eb 



cosy = 



(7) 



However, using semi realistic acceptance cuts of 
175° > > 5°, this method fails in the presence 
of finite width and radiative corrections, as 
can be seen from the fat histograms in figures 6 
and 7. 

On the other hand, it seems to be more promis- 
ing to try to reconstruct the neutrino momen- 
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Figure 6. Different methods of reconstructing the 
decay angle 0* in the presence of acceptance cuts 
(175° > 6» > 5°) and finite W± width. 
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Figure 7. Same as figure 6, but additionally in- 
cluding leading log initial state radiative correc- 
tions. 



turn from the missing transversal momentum, the 
mass shell condition = of the neutrino and 
the approximate mass shell condition of the lep- 
tonically decaying W: {p£ + p^e)^ = M^. If we 
add the further requirement that the squared sum 
of all reconstructed momenta should deviate no 
more than 10% from the squared CMS energy, 
we get the much improved hashed histograms in 
figures 6 and 7. They are in particular free from 
the characteristic shape distortions generated by 
the first method. 

5. Conclusions and outlook 

WOPPER is a fast, flexible and supported tool 
for physics at LEP-II and beyond. Compar- 
isons with other Monte Carlos [6] has shown that 
the approximations used in WOPPER (leading log 
radiation, resonant W^'s) can easily be controlled 
for experimentally relevant cuts. In the forthcom- 
ing releases, the following features will be added: 

• Anomalous couplings. This feature has 
been requested by experimentalists dur- 
ing this workshop again. Therefore this 
straightforward, if somewhat tedious, en- 
hancement will be installed. 

• Non electromagnetic radiative corrections. 
We will add the dominant contributions 
which go beyond the running QED cou- 
pling. The latter is of course already avail- 
able. 

• Coulomb singularity. We will add this fea- 
ture soon, again by popular demand^. For 
the time being, it will implemented in- 
clusively, without generation of the corre- 
sponding soft photons. 

• Improved photonic spectrum. Com- 
parison with complete 0{a) Monte Carlos 
shows that the LLA is an excellent approx- 
imation for inclusive distributions and lon- 
gitudinal spectra. There remains however a 
noticeable discrepancy at high photonic p'^, 

*At the time this contribution to the proceedings is writ- 
ten, the Coulomb singularity has already been imple- 
mented in WOPPER, Version 1.2. 
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which can not be described in LLA. This sit- 
uation will be improved, either by explicit 
inclusion of non leading terms or by phe- 
nomenological interpolation. 

• Final state radiation and interference 
terms. This will come somewhere further 
down the road. 

Due to the slowly varying nature of the e+e^ — » 
W~^W~ cross section, the forward branching al- 
gorithm [2] with hand crafted importance sam- 
pling as implemented in WOPPER vl.x is sufficient 
for most practical purposes. Wc shall however 
replace it by a backward branching algorithm, 
which will blend better with a general purpose 
Monte Carlo engine, where the cross sections can 
be replaced more easily [11]. 

6. Distribution policy 

WOPPER is distributed electronically over the In- 
ternet using the following channels: 

• The FORTRAN-77 sources (in PATCHY for- 
mat) can be obtained by anonymous Inter- 
net ftp from the host crunch, ikp.physik 
.th-darmstadt .de in the directories pub 
/ ohl/wopper/ old, pub/ ohl/wopper/pro 
and pub/ohl/wopper/new, corresponding 
to slightly outdated, current and experi- 
mental releases of WOPPER respectively. 

• The current status of WOPPER can be queried 
through the World Wide Web from the 
document http : //crunch . ikp . physik . th 
-darmstadt . de/monte-carlos . html. 

• Important announcements (new versions, 
fatal bugs, etc.) will be made through 
the mailing list wopper-cmnounceQcrunch . 
ikp . physik . th-darmstadt . de. Subscrip- 
tions should be mailed to wopper-announce 
-request at the same host. The purpose 
of wopper-announce is not general discus- 
sions of WOPPER, however, if there is inter- 
est among users, a companion list wopper 
-discuss can be created by the authors 
easily. 
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